Natural neighbor interpolation is a method of spatial interpolation, developed by Robin Sibson.[1] The method is based on Voronoi tessellation of a discrete set of spatial points. This has advantages over simpler methods of interpolation, such as nearest neighbor, in that it provides a more smooth approximation to the underlying "true" function.
The basic equation in 2D is:
where is the estimate at , are the weights and are the known data at . The natural neighbour method proposes a measure for the computation of the weights, and the selection of the interpolating neighbors
The natural neighbor method utilizes the change to the Voronoi tessellation to compute weights. The weights, , are by utilization of the area "stolen" from the surrounding points when inserting a new point into the tessellation. Each weight may be computed by dividing the section of the new tessellated region that lies within the tessellated region of each original neighboring tessellated polygon.